Kondo effect of impurity moments in <i-wave superconductors: 
Quantum phase transition and spectral properties 
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We discuss the dynamics of magnetic moments in d-wave superconductors, in particular we focus 
on moments induced by doping non-magnetic impurities into cuprates. The interaction of such 
moments with the Bogoliubov quasiparticles of the superconductors can be decribed by variants 
of the pseudogap Kondo model, characterized by a power-law density of states at the Fermi level. 
The Numerical Renormalization Group technique is employed to investigate this Kondo problem 
for realistic band structures and particle-hole asymmetries, both at zero and finite temperatures. 
In particular, we study the boundary quantum phase transition between the local-moment and 
the asymmetric strong-coupling phases, and argue that this transition has been observed in recent 
nuclear magnetic resonance experiments. We determine the spectral properties of both phases, the 
location of the critical point as function of Kondo coupling and doping, and discuss the quantum- 
critical cross-overs near this phase transition. In addition, the changes in the local density of states 
around the impurity are calculated as function of temperature, being relevant to scanning tunneling 
microscopy experiments. 



I. INTRODUCTION 

Impurities have proven to be a powerful probe for in- 
vestigating the bulk behavior of complex many-body sys- 
tems. In the field of high-temperature superconductivity, 
a variety of phenomena have been observed under dop- 
ing with magnetic as well as non-magnetic impurities: 
suppression of the superconducting critical temperature 
T c and an increase of the residual in-plane resistivityEra, 
damping of collective magnetic excitations (the so-called 
resonance modea), possible pinning of stripes and vor- 
tices, and impurity-related local bound, states seen in 
scanning tunneling microscopy (STM)qI3. 

A particularly interesting piece of physics is the mag- 
netism of impurities which are substituted for Cu ions. 
Experiments have been performed with magnetic spin- 
1 (Ni) as well as non-magnetic spin-0 (Zn, Li) impuri- 
ties. Whereas the spin-1 impurities naturally carry an 
on-site moment which is expected to give rise to some 
kind of Kondo physics, the behavior of non-magnetic im- 
purities is more surprising. A series of hppjtiful nuclear 
magnetic resonance (NMR) experimentsErE-3 have clearly 
shown that each impurity, despite having no on-site spin, 
induces a local 5=1/2 moment on the neighboring Cu 
ions at intermediate energy scales. In particular, the local 
susceptibility associated with Li ions in YBa2Cu306+a; 
(YBCO) has been found to show Curie- Weiss-like behav- 
ior, with a strongly doping dependent Weiss temperature 
which appears to vanish for strongly underdoped sam- 
ples - this implies the existence of free moments in the 
underdoped regime down to temperatures of order 1 K. 
Microscopic perspectives on the formation of these local 
impurity moments will be discussed in the next section. 

Accepting that non-magnetic impurities induce local 
moments, these moments are expected to interact with 



the elementary excitations of the bulk materialS In this 
paper, we will be mainly concerned with a d-wave super- 
conducting bulk state, i.e., with temperatures T below 
the superconducting T c : the relevant low-energy excita- 
tions of the d-wave superconductor are fermionic Bogoli- 
ubov quasiparticles and bosonic antiferromagnetic spin 
fluctuations (as seen in neutron scattering); quantum 
phase fluctuations which may become important near a 
superconductor-insulator quantum phase transition will 
not be considered here. The interactions of impurity mo- 
ments with fermionic and bosonic degrees of freedom can 
be discussed separately, as they lead to distinct phenom- 
ena associated with fermionic and bosonic Kondo models. 
Also, the energy scales of these two phenomena appear 
well separated, as the fermionic Kondo temperature in 
the superconducting state is below 100 K whereas the 
spin fluctuation energy scale is given by the energy of 
the "resonance mode"H being 40 meV at optimal doping. 

The interaction of an impurity spin with bosonic spin- 
1 fluctuations of a quantum disordered aniifcppmagnet 
is described by a "bosonic Kondo modeF'Ejta. It has 
been studied recently in the context of impurities in 
two-dimensional nearly-critical antiferromagnets, where 
it leads to a new (2+l)-dimensional boundary quan- 
tum field .theory with a number of interesting universal 
properties!!!! If the bulk state has a finite spin gap then 
the impurity moment is not screened, i.e., it contributes 
with a Curie term to the impurity susceptibility. How- 
ever, the moment is spatially "smeared" over a length 
scale given by the magnetic correlation length, i.e., part 
of the moment is carried by neighboring spins. Another 
result of Ref. 16 was that a finite concentration of impu- 
rity moments leads to universal damping of the collective 
spin-1 mode of the bulk magnet. Together with recent 
neutron scattering experimental this provides further ev- 
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idence for 5 = 1/2 induced local moments near Zti sites 
in cuprate superconductors: it has been arguedtj that 
these moments are required to explain the strong effects 
of a small concentration of Zn impurities on the "reso- 
nance peak" in the spin dynamic structure factor. 

In this paper, we will focus on the interaction of the 
impurity moment with the fermionic bulk degrees of 
freedom, namely the quasiparticles of the <i-wave su- 
perconductor. Importantly, the single-particle density 
of states (DOS) of the Bogoliubov quasiparticles van- 
ishes at the Fermi energy, and so we expect a number 
of features quite distinct from the usual Kondo effect in 
metalsliia. The Kondo effect in systems where the host 
DOS follows a power-law near the Fermi level, p(e) ~ |e| r 
(r > 0), has been studied quite ext^si«cly in the con- 
text of-thc. jpspudogap Kondo mode]E3~L3. A number of 
studiesEjOSt3 including the initial work by Withoff and 
Fradkin employed a slave-boson large- N techniqueEil; fur- 
ther progress and insight came fromjuimcrical rcnormal- 
ization group (NRG) calculationsElLJ, the lpc.al moment 
approachEZI, and a dynamic large- N methodEl The gen- 
eral picture arising from these studies is that there ex- 
ists a zero-temperature (boundary) phase transition at 
a critical Kondo coupling, J c , below which the impurity 
spin is unscreened even at lowest temperatures. Also, the 
behavior depends sensitively on the presence or absence 
of particle-hole symmetry: in the particle-hole symmet- 
ric case there is no complete screening even for Kondo 
couplings J > J c . A comprehensive discussion of possi- 
ble fixed points and their thermodynamic properties has 
been given by Gonzalez-Buxton and IngersentES based 
on the NRG approach. Fig. |l| summarizes their findings 
for r > 1/2, including the case r = 1 relevant for d-wave 
superconductors (see also Fig. 16 of Ref. p6|) . 

Spectral properties of the pseudogap Kondo model (for 
the relevant case of SU(2) spin symmetry) have so far 
been-ipvestigated for the particle-hole symmetric case 
onlycjcZL However, interesting properties are to be ex- 
pected especially fo^the-asymmetric ca f%, considering 
the results of statidfJ'Eirc3 and dynamicEa large- N ap- 
proaches. 

The study of spectral properties of Kondo impurities 
in metals as well as in superconductors has been revived 
due to recent advances in scanning tunneling microscopy 
(STM). Tunneling to clean surfaces of materials with 
bulk or surface impurites has made it possible to directly 
observe local spectral properties associated with impuri- 
ties. In particular, Kondo impurities on metal surfaces 
have been found to give rise-.to, a dip or Fano-like line- 
shape in the tunneling DOSE3H. In contrast, in cuprate 
d-wave superconductors peak structures in the local DOS 
close to the Fermi energy have been found for both 
Zn and Ni impurities U~U doped into Bi 2 Sr 2 CaCu208 + 5 
(BSCCO). Whereas the results for Ni appear to be well 
described by a model with combined potential and static 
magnetic (spin-dependent) scatteringcj, the signal of Zn 
impurities - a huge tunneling peak at energies of 1-2 meV 
- is more puzzling. It is tempting to identify this peak 
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FIG. 1. Schematic zero-temperature renormaliza- 

tion-group flow diagramEj for the pseudogap Kondo model 
with power law density of states p(e) ~ |e| r with r > 1/2. 
The horizontal axis denotes the Kondo coupling, the vertical 
axis denotes particle-hole asymmetry (parametrized, e.g., by 
the potential scattering strength U for a particle-hole sym- 
metric host band structure). The model has two stable fixed 
points (solid dots): the weak-coupling local moment (LM) 
fixed point where the impurity is essentially decoupled from 
the band, and the asymmetric strong coupling (ASC) fixed 
point where the impurity moment is fully screened. The thick 
line represents the second-order phase transition between the 
two zero-temperature phases, the open dot is the critical fixed 
point. Note that there is no screening in the particle-hole 
symmetric case even for infinitely strong Kondo coupling. 

with a quasiJaound state in a purely potential scatter- 
ing modelEa~E3, but such a state appears at low energies 
only for a range of very larg e . poten tial values depend- 
ing upon microscopic detail£j'E3El Furthermore, the 
spatial dependence of the zero bias peak is surprising - 
further discussion in Sec. [v| - and the observed spatially 
integrated spectrum is asymmetric between positive and 
negative biasu, whilp-,the potential model predicts apt. 
proximate symmetrycll. Recently, it has been proposed^ 
that some of the properties of the Zn resonance can be ex- 
plained by taking into account the Kondo spin dynamics 
of the magnetic moment induced by Zn. The screening of 
this moment by the Boboliubov quasiparticles provides 
a natural low-energy scale (the Kondo temperature) ex- 
plaining the energetic location of the peak seen in STM. 

The purpose of this paper is twofold: On one hand we 
want to study the fermionic pseudogap Kondo effect us- 
ing realistic band structures for the cuprates. We shall 
show that the transition between a free moment and a 
screened moment in the pseudogap Kondo model can be 
driven by varying doping in the high-T c compounds - 
this provides an explanation for the strong doping de- 
pendence .of. the Kondo (Weiss) temperature observed 
in NMROEj. On the other hand we want to examine 
a number of properties of the pseudogap Kondo model 
close to this transition and relate them to experimental 
findings, in particular STM measurements. 

The remaining part of the paper is organized as fol- 
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lows. In Sec. || we briefly discuss the issue of properly 
modelling a non-magnetic impurity in cuprates - we will 
motivate the model consisting of a Kondo and a potential 
scattering term, a nd fu rther discuss a number of theoret- 
ical aspects in Sec. III. In Sec. IV we turn to the simplest 
model, namely a point-like magnetic impurity with a re- 
alistic host density of states (taking dispersion and gap 
data of actual cuprate materials). Using the numerical 
renormalization group method, we will discuss static and 
dynamic properties of the two stable fixed points as well 
as of the quantum critical regime. We will determine 
the critical coupling constant and the cross-over scale as 
function of doping, and make contact with the transi- 
tion seen in NMR experiments. In Sec. [v| we extend the 
impurity modelling to a spatially distributed magnetic 
moment and additional potential scattering and present 
in particular the resulting local conduction electron DOS 
which is relevant for STM experiments. A brief discus- 
sion of open issues will close the paper. 

Parts of the discussion in Sec. || have already been 
given in Refs. 42 and 43, but we decided to include them 
here to keep the paper self-contained. Readers mainly 
interested in the experimental consequences for cuprates 



might skip t he technical issues in Sec. Ill, and have a l ook 
at Sec. IV C (especially Fig. |qJ> , as well as Sees. IV B 2 and 
which are related to STM observations. 



II. EFFECTIVE MODELS FOR NON-MAGNETIC 
IMPURITIES 

As NMR experiments show, non-magnetic impurities 
can induce local moments in correlated hosts. In this 
section we will briefly discuss theoretical perspectives on 
this phenomenon. 

To describe a non-magnetic impurity in a strongly cor- 
related superconducting system, two approaches appear 
possible, (i) One starts with a basic model for the bulk 
system containing correlation terms, e.g., a one-band 
or three-band Hubbard model with strong on-site inter- 
action, and describes the non-magnetic impurity by a 
potential scattering term as suggested by its chemistry. 
This approach is notoriously difficult, as a strongly cor- 
related model (which is hard to deal with even without 
impurity) has to be treated including the impurity and in 
the low-temperature limit. This task is almost impossible 
for current numerical .techniques such as Monte-Carlo or 
exact diagonalizationa, and analytical approaches suf- 
fer from uncontrolled approximations, (ii) One separates 
the questions of moment formation and interaction of the 
moment with the bulk system. To describe the latter, 
one can employ effective models for both the impurity 
and the bulk system. The model for the bulk system 
has to contain only the minimal ingredients to describe 
the desired low-temperature physics - in the present case 
this is a d-wave BCS model which does not contain corre- 
lation effects other than the Hartree-Fock pairing term. 



The impurity term of the effective model has to account 
for the correlation effects related to the introduction of 
the impurity - in the case of interest it has (at least) to 
contain the magnetic moment and a scattering potential. 

In this discussion and the rest of the paper, we will 
adopt the approach (ii), see also Ref. Therefore, we 
will not attempt a rigorous derivation showing the ex- 
istence of the magnetic moment near the non-magnetic 
impurity, but instead give some plausible arguments for 
its appearance. If we consider an undoped paramagnetic 
Mott insulator as reference system at T — 0, the moment, 
formation is easily seen to arise from breaking singletg3'c3 
by removing spins. These unpaired spins will remain in 
the vicinity of the impurities if the host antifcrromagnet 
has confined spinons, i.e., elementary spin-1 excitations. 
In this picture each impurity can be understood as lo- 
calized "holon" which binds the moment of a S = 1/2 
"spinon" . By continuity, the described mechanism is ex- 
pected to be effective also at small, but finite hole dop- 
ing. At larger doping, a related picture can be developed 
by analogy with the theorycj for moment formation in 
the disordered metallic state of Si:P — small variations in 
the potential combine with strong local correlations to 
induce very localised spin excitations. Other theoreti- 
cal perspectivesc3~E3 on local moment formation in the 
cuprates have also been given - most of them are based on 
the related idea that strong potential scatterers can form 
quasi-bound states at the impurity sites near the Fermi 
levella. Accounting for the strong local Coulomb repul- 
sion, each bound state will capture only a single electron, 
and the low energy physics will again be described by an 
impurity spin model. 

It is evident that the magnetic properties associated 
with the induced moment will strongly depend on the 
hole doping level. The undoped limit, i.e., the undoped 
paramagnetic Mott insulator, has a free 5 = 1/2 moment 
near each impurity, while the strongly doped limit (where 
electronic correlations are presumably weak) does not. 
From continuity one expects that the free T = moment 
will survive in the superconducting state for a finite range 
of doping; spin quantization suggests that the size of the 
moment as measured by the Curie term in the suscepti- 
bility is always S — 1/2 independent of doping, i.e., each 
moment will contribute a divergent Curie susceptibility 
~ 1/4T even in the (underdoped) superconducting state. 
Then, a quantum critical point separates the weak and 
strong doping limits. On the strong doping side of this 
quantum critical point, the moment is Kondo screened 
as T — ► 0, and such a regime is continuously connected 
to a regime, at higher doping, where the moment does 
not even form at intermediate T. We wish to argue that 
this quantum phase transition is precisely the transition 
present in the pseudogap Kondo model, and this claim is 



supported by quantitative calculations in Sec. IV 



We note that the pictures of the moment formation 
quoted above do not give reliable information about 
the parameters of a possible effective impurity model. 
NMR suggests that the magnetic moment is spatially 
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distributed among the Cu sites near the impurity (but 
fluctuates as a single entity). However, the precise mi- 
croscopic form of the coupling between the spin moment 
and the conduction electrons is not known - one has to 
keep in mind that the moment is formed by a particu- 
lar (bound) state of the conduction electrons near the 
impurity, and it interacts with the other states (i.e. lin- 
ear combinations) of the same conduction electrons. In 
absence of a precise knowledge of the appropriate micro- 
scopic model we will employ some simple Anderson- and 
Kondo-like models for the impurities which are discussed 
in the following section. 

III. MICROSCOPIC MODELS AND METHODS 

For the remainder of this paper we restrict our atten- 
tion to effective low-energy models consisting of a BCS- 
like d-wave superconductor, a scattering potential, and a 
magnetic moment, H = Hbcs + Wpot + W mag . In partic- 
ular, we neglect the strong short-range magnetic corre- 
lations of the bulk material. As mentioned in the intro- 
duction, these low-energy magnetic fluctuations couple to 
the impurity moment leading to a bosonic Kondo model. 
Due to the non-zero spin gap, the main effect of this cou- 
pling is an additional spatial "smearing" of each impurity 
moment, which can be absorbed in the effective impurity 
model. 

We assume that the cuprate superconductor can be 
described within BCS theory, so we use 

W B cs = X>k[(£k-M)r z +A k T 3 ']v]> k . (1) 

k 

Here vf^ = (ck|! c Lkj.) ^ s a Nambu spinor at momentum 
k = (k x , k y ) (cka annihilates an electron with spin a on a 
Cu 3d orbital), t x ' v ' z are Pauli matrices in particle-hole 
space, and fj, is the chemical potential. For the kinetic en- 
ergy, £k, we will use a tight- binding form which includes 
nearest neighbor hopping as well as longer-range hopping 
processesEU, while we assume a d-wave form for the BCS 
pairing function Ak = (A /2)(cos k x — cosk y ). 

In the numerical calculations of this paper we will ne- 
glect order parameter relaxation due to the impurity. 
This effect might play a role for strong scattering po- 
tentials, in that it suppresses the magnitude of the gap 
and leads to a finite zero-energy DOS. However, STM 
tunneling experimentstffl show that the local change in 
the gap size is rather small. Furthermore, the induced 
zero-energy DOS will also be small, and the resulting (fi- 
nite) Kondo temperature is exponentially suppressed. So 
we expect results including order parameter relaxation to 
be similar to ours below, at least regarding the magnetic 
properties of the impurity in the experimentally accessi- 
ble temperature range. 



A. Anderson vs. Kondo model 

An effective model for a (magnetic or non-magnetic) 
impurity consists of a potential scattering term and a 
magnetic term. The potential scattering term is assumed 
to be localized at the impurity site: 

n P ot = uJ2^o a - (2) 

a 

The problem of a single scatterer can be solved exactly, 
with the standard result for the 'J Green's function in 
Nambu notation: 

G(r, r', w) = G°(r - r', «) - UG°(r - r , cu n ) 

xt z [1 + [^((O.O^u/ItT^Vo - r»; (3) 

G° is the Green's function of the host N s G a (r,oj) = 
£ k e lk -> - (e k - fi)r z - A k T x }~\ and r = (0,0) is 
the scattering site. 

The magnetic moment can be either coupled to a single 
site (point-like impurity) or be spatially distributed. For 
the point-like magnetic impurity coupled to a single site 
r = (0, 0) of a (metallic or superconducting) host the 
most general form is a single-impurity Anderson model: 

ftmag = ^2 V 0(ft C 0<r + h.C.) 

• 'fY^LJ^ ■ '<//;•»/.• ( 4 ) 

cr 

In the so-called Kondo limit, Vq — > oo, ey — > — oo, U — * 
oo, charge fluctuations on the impurity orbital are frozen, 
and the Anderson model can be mapped onto a Kondo 
model (Schricffer- Wolff transformation): 

Hmag = • S (5) 

where s = N^ 1 X)kk' Q/ 3 c ka5°W3 c k'/5 is the conduction 
band spin operator at the impurity site i"o = (0, 0), N the 
number of lattice sites, and J = 2V r 2 (l/|e/| + l/|t r /+e/|). 
For Uf/2 =/= —e/ this mapping also introduces an addi- 
tional potential scattering term which can be absorbed 
in U (§. 

Both models allow for a straightforward generalization 
to spatially distributed impurities, which are, however, 
no longer equivalent! The Anderson model 

Hmag = ^ VniflcRa + h.C.) 

+ e/^/]/a + C//n /T n /i (6) 

cr 

is easily seen to describe an Anderson impurity coupled 
to a single linear combination of conduction electrons 
on the sites R, c = J2r ^R c r/^ where V 2 = J2r ^r- 
In the Kondo limit, such a model maps onto a Kondo 
model with non-local Kondo couplings. In contrast, the 
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straightforward generalization of the Kondo model ([3|) to 
an extended impurity reads 

Hmag = ^ JrS ■ Sr , (7) 
R 

this model represents a multi-channel Kondo modelElEl. 
The impurity spin couples to all possible linear combina- 
tions of the conduction electrons on the R sites. In other 
words, a spatially extended Kondo impurity generically 
couples to different angular momentum channels of the 
conduction electrons. 

The relation between the models (||) and (Q) has been 
recently discussed in Ref. p4^_ It has been argued that 
in an Anderson model like (g) the coupling to a corre- 
lated host opens new screening channels, and the effec- 
tive model will be a multi-channel Kondo model (Q) . The 
main idea is that conduction band correlations reduce 
charge fluctuations in the host, and if an electron hops 
onto the impurity site, it has to hop back preferentially 
onto the same conduction band site where it came from. 

We note that the different screening channels gener- 
ated by such a mechanism are certainly not equivalent, 
and the low-T behavior will be determined by a single- 
channel fixed point. (However, this does not exclude that 
genuine multi-channel physics may be realized at inter- 
mediate temperatures.) We will employ a Kondo model 
of the form ((?]) in Sec. when discussing STM experi- 
ments on Zn in BSCCO. 



B. Influence of superconducting correlations 

A Hamiltonian of the form Hbcs +H m ag does not only 
describe a host with a pseudogap density of states, it also 
contains superconducting correlations. From diagram- 
matic pertubation theory it is clear that the Kondo prob- 
lem in a normal-state system with a power- law DOS (e.g. 
a semi-metal) is not equivalent to the Kondo problem in 
a superconductor with the same DOS - in other words, 
anomalous Green's function contribute to the screening. 

For the special case of a point-like Kondo impurity io. 
an s-wave superconductor, however, it has been shownO 
that the problem with anomalous Green's functions can 
be mapped onto a Kondo problem in a normal-state 
system with modified gapped DOS, containing an ad- 
ditional particle- hole asymmetry (i.e., a potential scat- 
tering term), but without anomalous Green's functions. 

In d-wave superconductors such a mapping has not 
been achieved so far, and the situation is less clear be- 
cause the anomalous Green's functions of the d-wave su- 
perconductor provide a coupling between different an- 
gular momentum channels. Therefore, it is expected 
that the Kondo problem in a d-wave superconductor in- 
volves infinitely many bands (angular momenta) even for 
a point-like impurity. 

Various approximate methods which have been em- 
ployed for solving the Kondo problem in a d-wave super- 



conductor handle the anomalous Green's functions dif- 
ferently. Most calculations turn out to be insensitive to 
superconducting correlations in the host material, i.e., 
they treat a semi-metal rather than a superconductor. 

For instance, in the standard slave-boson approach, 
anomalous terms are formally included, however, for a 
point-like impurity in a <i-wave superconductor their con- 
tributions are easily seen to drop out. The standard 
implementation of the numerical renormalization group 
method uses the host density of states as the only input, 
and anomalous terms are ignored. Based on the results 
for s-wave superconductors, we nevertheless expect that 
the results obtained by NRG are at least qualitatively 
correct in a generic particle-hole asymmetric situation, 
as the main influence of the anomalous Green's functions 
is possibly an additional particle-hole asymmetry. 

C. Numerical renormalization group method 

The NRG method has been developed by Wilson for 
the investigation of the Kondo modeled. Due to the log- 
arithmic discretization of the conduction band, it is able 
to access arbitrarily low energy scales which is essential 
also for the model discussed in this paper. (However, due 
to the logarithmic discretization it is not possible to re- 
solve sharp structures at high energies.) The NRG is a 
non-perturbative method as the impurity site including 
the strong local Coulomb interaction is treated exactly. 

The NRG method has been generalized to a number of 
different impurity models^such as the pseudogap Kondo 
and Anderson modelsEJo. For these models, the main 
modification is the mapping of a non-constant conduction 
electron density of states onto a semi-infinite chain form 
(see, e.g., Ref. p4| ). 

Recently the NRG method has been extended to the 
calculation of dynamic properties at finite temperatures, 
we refer the reader to Ref. |5?] for details. We should 
mention here that information on spectral quantities for 
frequencies lo <C T cannot be determined from the NRG 
method; the resolution is limited by the temperature and 
possible structures at frequencies much smaller than the 
temperature are missed by this approach. However, this 
causes no problem when comparing with, e.g., STM ex- 
periments, because the finite temperature in both the 
sample and the tip introduces a thermal broadening, such 
that structures on scales smaller than T will not be re- 
solved. 

The NRG calculations presented in this paper are 
mostly performed with a discretization parameter A = 2 
and 600 or 800 levels. 



D. Slave-boson mean-field approximation 

The slave-boson mean-field approximation!^] to the 
Kondo impurity model is a large- A approach, i.e., it be- 
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comes exact in the limit of large spin degeneracy N — > oo. 
It is known that this method is not quantitatively accu- 
rate for the physical case N = 2, and has numerous ar- 
tifacts, both at finite temperatures in general as well as 
near the quantum-critical point in the pseudogap Kondo 
model. However, it is expected to capture the qualitative 
physics in a Kondo screened phase. 

The method has been applied ia a-aumber of papers to 
the pseudogap Kondo modelc2rEilEj'E2l, and we will men- 
tion some results below. 



IV. POINT-LIKE MAGNETIC IMPURITIES 

In this section, we discuss general properties of the 
asymmetric pseudogap Kondo model, in particular the 
quantum phase transition between the local-moment and 
the asymmetric-strong-coupling phase. Since we are 
mainly interested in the universal properties near this 
transition, we restrict our attention here to a S = \ im- 
purity coupled to a pseudogap band structure, and give 
most results for general values of the pseudogap exponent 
1/2 < r < 1 (some results are valid also for < r < 1/2). 
We will obtain numerical results using NRG, where the 
DOS is the only host quantity entering the calculation. 
For technical purposes, namely a convenient calculation 
of the impurity spectral function (or T matrix) , the impu- 
rity is modelled by an Anderson Hamiltonian (Q) (instead 
of a Kondo Hamiltonian) with Uf = — 2ey, and param- 
eters \Uf\, \ef\ ^-.bandwidth which places the model in 
the Kondo limitES. 

When specifying the results for the high-T c com- 
pounds, we consider an impurity coupled locally to 
a single site of the host <i-wave superconductor. To 
obtain quantitative estimates for characteristic ener- 
gies, we use tight-binding hopping parameters ex- 
tracted from normal-state angle-resolved photoemission 
(ARPES) datatil and sizes of the <i-wave gap as obtained 
from a number of tunneling, .and ARPES experiments in 
the superconducting stateE3E3. 



A. Zero-temperature phase transition and static 
properties 

As_dis£ussed in a number of papers over the last 
yearaHTEj the pseudogap Kondo model (with power-law 
DOS p(e) ~ |e| r ) shows a phase transition as a function 
of the Kondo coupling. For pseudogap exponents r > 1/2 
there are two stable zero-temperature fixed pointsc3, 
namely (i) a local moment (LM) regime is reached for 
weak initial coupling - here J flows to zero, and the im- 
purity is unscreened - and (ii) an asymmetric strong- 
coupling (ASC) regime which is reached only for suffi- 
ciently large Kondo coupling and particle-hole symmetry 
breaking. The renormalization group flow is sketched in 
Fig. | 



Particle-hole asymmetry is a relevant parameter in the 
pseudogap Kondo problem, as seen in Fig. ||. For a 
particle-hole symmetric host band structure, it can be 
parametrized by the amount of potential scattering at 
the impurity. In the general case of an asymmetric host, 
particle-hole asymmetry can still be cast into a single 
number, in the case of a point-like Kondo impurity the 
real part of the local host Green's function (after inclu- 
sion of potential scattering) is a suitable choice. We will 
see below that, e.g., the location of the "Kondo" peak 
will depend on the sign of this overall particle-hole asym- 
metry. Concerning the cuprates, this overall sign cannot 
be extracted from experiments or theory, as it depends 
on details of the two-dimensional band structure which 
place the van-Hove singularity either below or above the 
Fermi level. (Remarkably, this van-Hove singularity has 
never been observed experimentally.) 

For the following discussion, it is useful to define sus- 
ceptibilities describing the response to external magnetic 
fields. We allow for a space-dependent field H u coupled 
to the host and a local field H lmv at the impurity, 



H 



mag 



H 



mag 



(8) 



With these definitions, a space- independent uniform field 
applied to the whole system corresponds to H u (i) = 
= H. Response functions can be defined from 



derivatives of the free energy, F 
as follows: 



-TbxZ (k B = 1) 



Xu,u(*; O — 
Xu,imp(*) — 
Ximp,imp = 



T 



S 2 InZ 



3 5H ua (i)SH ua {i>) 
T S 2 \nZ 
3 SH ua (i)SHi mp ^ Q 
T S 2 lnZ 

3 fiffimp ^imp,a 



(9) 



From these quantities we can define various observables, 
starting with the impurity contribution to the total sus- 
ceptibility, 



Ximp(T) = x 



imp, imp 



2 ^ ^ Xu,imp(^) 



EMM')-^ k (v')) 



(10) 



where Xul| k i s the susceptibility of the bulk system in 
absence of impurities. The local impurity susceptibility, 
i.e., the response of the impurity spin to a local field, is 
given by 



XiocCO = x 



imp, imp j 



(ii) 



which is equivalent to the zero-frequency impurity spin 
autocorrelation function. Importantly, NMR Knight 



G 



shift experiments measure a different quantity - there 
a uniform field is applied to the whole sample, and the 
local response is given by 



Xnmr(T) = Xi 



imp, imp 



imp 



(i). 



(12) 



Note that sometimes this quantity is referred to as "local 
susceptibility", e.g., in Ref. ||^. 

After having defined the susceptibilities, we return to 
the pseudogap Kondo problem. We briefly recall the im- 
purity thermodynamic properties of the relevant stable 
fixed points (LM and ASC) - these have been given, 
e.g., in Ref. Eq. The local-moment fixed point has all 
characteristics of a decoupled spin- 1/2: susceptibility 
^Ximp = 1/4, entropy 5i mp = In 2, and specific heat 
Cimp = 0. At the asymmetric strong-coupling fixed point 
the impurity spin is fully quenched: T^imp = 0, S\ mp — 1 ' 



and Ci 



imp 



0. 



Remarkably, the leading correctionsEa 
here are A(T Xlmp ), A5 imp , AC imp cx T 2r for < r < 1, 
in other words, the impurity susceptibility Ximp(T') van- 
ishes as T — > for r > 1/2. 

We note that there exists a third stable fixed point of 
the pseudogap Kondo model, namely a symmetric strong 
coupling fixed point. It is reached for pseudogap expo- 
nents r < 1/2 in the particle- hole symmetric case for 
sufficiently strong Kondo coupling; its thermodynamic 
properties are-spmewhat different from the ASC behavior 
quoted abovecn However, the real materials are gener- 
ically particle-hole asymmetric, and the physics of the 
particle-hole symmetric strong-coupling fixed point is of 
no relevance here. 



B. Quantum-critical and cross-over behavior 

In this section, we focus on the properties of the impu- 
rity model in the vicinity of the quantum-critical point 
separating the local moment from the asymmetric strong 
coupling phase. (Critical properties for the particle-hole 
symmetric pseudogap model are different, and have been 
examined, e.g., in Refs. ^ and ^3 .) 

In the vicinity of the critical point, one can define an 
energy scale T* , which vanishes at the transition, and de- 
fines the cross-over energy above which quantum-critical 
behavior is observedO. If T* is much smaller than both 
the bandwidth and the superconducting gap we expect all 
observables to show scaling behavior as function oiT/T* 
and w/T*. Importantly, for the special case of r = 1, log- 
arithmic corrections occur - strictly speaking, these cor- 
rections invalidate one-parameter scaling behavior since 
a second (high-energy) scale A (of order bandwidth) en- 
ters in the form of ln(T/A) or ln(w/A). However, for 
most practical purposes the logarithmic corrections are 
small, and in any case do not change the leading power 
law behavior. We have extracted the logarithmic cor- 
rections from the NRG results for some observables, and 
they are quoted below. 
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FIG. 2. Cross-over diagram of the pseudogap Kondo 
model as function of the Kondo coupling J and temperature 
T, calculated with a point- like impurity and a host band struc- 
ture corresponding to a cuprate superconductor at optimal 
doping. The dot on the horizontal axis denotes the quantum 
phase transition between LM and ASC phases, the dashed 
lines are cross-overs, indicating T*. Both quantum-critical 
and local-moment regimes show a impurity susceptibility di- 
verging as 1/T, whereas in the asymmetric strong coupling 
region the impurity susceptibility vanishes oc T. As we will 
see in Sec. [V], the critical Kondo coupling becomes signifi- 
cantly smaller for a spatially distributed impurity moment. 



We start with the dependence of the cross-over scale 
T* on the reduced coupling, j = (J — J c )/J c , measuring 
the distance from the critical point at J c . Our numerical 
results obtained by NRG are consistent with 



l/r 



(13) 



for 1/2 < r < 1 on both sides of the transition. In Fig. g 
we show numerical results for T* for a high-T c compound 
band structure at optimal doping. In Fig. || it is clearly 
seen that even for Kondo couplings well away from J c , 
the cross-over scale T* is small (20 K), therefore impurity 
quantum-critical properties may easily be observed in the 
cuprates in a rather large temperature range without fine 
tuning. Parenthetically, we note that the cross-over scale 
T* is in general defined only up to a prefactor of order 
unity; in Fig. || the numerical value of T* was taken from 
the location of the maximum in the T matrix spectral 
density, see below. 

For r = 1 the cross-over temperature T* vanishes lin- 
early with the distance to the critical point (up to log- 
arithmic corrections). ._ This is in agreement with ear- 
lier NRG calculations^] as well as with the dynamic 
large-A" approacliEx However, it is at variance to 
the results obtained within the slave-boson mean-field 
approximationEI which predicts an essential singularity 
for r = 1, i.e., an exponentially vanishing scale T* near 
the transition. We believe that the NRG prediction of a 
linearly vanishing T* is correct (as the NRG is reliable in 
capturing asymptotic low-energy physics), and that the 
exponentially vanishing scale is an artifact of the slave- 
boson method. 
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FIG. 3. a) NRG results for the temperature dependence 
of TXimp, for r = 1, using a band structure corresponding to 
optimally doped cuprates. The different curves correspond to 
different values of the Kondo coupling J which are very close 
to the critical coupling J c « O.lleV. b) Scaling function 
<E>Xim P (T/r*) extracted from NRG results, for r = 1. (Pos- 
sible logarithmic corrections violating scaling are very small 
here.) The value of T* has been taken fro m the peak position 
in the impurity spectral density, see Sec. [VB2. 



1. Susceptibility 

As a first observable we consider the impurity suscep- 
tibility XimpfT)- As a response function associated with 
a conserved, quantity, it cannot acquire an anomalous 
dimensionEj, therefore we expect precisely at the criti- 
cal point a temperature dependence as 



Ximp(r, J=J C ) 



c 

T 



(14) 



Here C is a universal number, depending only on r. In- 
terestingly, the behavior ( |l4| ) can be interpreted as the 
Curie response of an irrational spin. (C = 1/4 would cor- 
respond to the usual Curie response of a free spin 1/2.) 
Numerical data for Tximp are shown in Fig. |^. In par- 
ticular, close to J c , there is a plateau in Tximp, which 
allows to extract the universal constant C(r = 1) w 0.14 
in (|l4|). This number is somewhat smaller than the value 
0.164 quoted by Gonzalez-Buxton and Ingersent (Fig. 14 
of Ref. |26|), this might be due to discretization or cut-off 
effects in the numerics (we have not attempted a careful 
finite-size analysis here). 

The local susceptibility (|ll]) is not protected by a con- 
servation law, and can acquire an anomalous exponent 



?7i oc . This implies a power law 

Xioc(T, J = J C ) (xT" 1 "- 1 



(15) 



at the critical point. 

For both finite T and finite ( J — J c ), the low-energy be- 
havior is completely determined by the cross-over energy 
scale T* and the temperature T itself (up to logarithmic 
corrections for r = 1). Then, the susceptibilities can be 
described by universal cross-over functions, 



Ximp(T) = ^® Ximp (J^) 



T 2 



(16) 



with <I>Ximp(» = ^Xim P {x) / x 2r , and 



B 



T 



J 1 !— 'Jloc Xloc \ r £* 



T 



(17) 



with $X\oc{ x ) = 2;' ?loc< i ) Xioc(a ; ); and B an amplitude pref- 
actor. The equations (|l6| ) and ( |l7| ) define universal scal- 
ing functions $ of the reduced temperature T/T*; the $ 
are, of course, different for different r and for both sides 
of the quantum phase transition, i.e., we have to distin- 
guish $ LM and $ ASC . The scaling functions $Xim P for 
r = 1 are shown in Fig. |^b - here our present numerics 
is not accurate enough to observe possible logarithmic 
corrections to scaling. Also, a reliable fitting of Ximp to 
a Curie- Weiss law for a temperature range above T* re- 
quires a finer discretization in the NRG procedure, and 
has not been attempted here. It is clear that the Weiss 
temperature will be approximately given by the cross- 
over temperature T* (defined by the T matrix peak), 
but we cannot give a reliable estimate for the ratio of the 
two. 

We briefly discuss the asymptotics of the scaling func- 
tions $ x . Spin quantization requires $^ p (0) = 1/4 and 

$ Xilp(°) = °> tncrcforc the zero-temperature limit of 
Tximp is fixed to 1/4 throughout the whole local-moment 
phase (and to in the whole ASC phase). The uni- 
versal Curie response at the critical point Jl4| ) implies 

$ xiip(°°) = $ x^ P (°°) = c - From Xim P (r -> °) K t2,_1 

in the ASC phase we deduce $^ c (a;) oc x 2r for small 
x. This immediately gives Xim P {T = 0) oc T 2 ^ 1 /T* 2r oc 
T 2r ^ 1 /(J — J c ) 2 for J > J c . However, the real super- 
conductor will have a small residual DOS at the im- 
purity site due to order parameter relaxation, this will 
will lead to a finite impurity susceptibility in ther-aero- 
temperature limit for J > J c . Finally, NRG datacS in- 
dicate that xioc(T — > 0) oc 1/T in the LM regime; this 
gives us (a;) oc x~ Vl ° c , and therefore Tx\ oc (T — > 0) oc 



T* ni °" oc {J-J c ) 



Vloc/r 



in the LM phase. Our preliminary 
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NRG calculations of Tx\ oc allowed to extract the value 
of 7710c, we found 77i oc w 0.05 for r = 1. 

At this point a few remarks about experimental sus- 
ceptibility measurements and the possibility to observe 
the described quantum-critical behavior are in order. As 
detailed above, the only T dependence with a non-trivial 
power law can be expected in the local susceptibility 
X\oc{T)- However, a direct measurement of xi oc requires 
a local field, and is perhaps only possible using muon spin 
resonance. In contrast, commonly used NMR techniques 
probe Xnmr ( |l2| ) , and this quantity shows either a Curie 
law or saturates to a constant at low T. Measurements of 
the total susceptibility in the local-moment regime should 
observe Curie behavior with the "full" prefactor corre- 
sponding to spin-1/2 per impurity, i.e., Ximp = Mmp/4T, 
where iV; mp is the number of impurities. Such measure- 
ments have been performed on Zn-doped YBCOE3, and 
the results indicate a somewhat smaller value of the mo- 
ment per impurity, which on the other hand depends 
only weakly on doping. Two explanations for the de- 
viation from the S — 1/2 moment appear to be possi- 
ble: Either there are subtle cancellation effects involved 
in the SQUID susceptibility measurements (one also has 
to keep in mind that the subtraction of the bulk sus- 
ceptibility leads to rather large uncertainties in Ximp), 
or the observed reduction of the moment reflects the 
quantum-critical behavior (|l4|). In the latter case, one 
should expect a cross-over between the Curie law (|1J) to 
a Curie law with the "full" prefactor at low enough tem- 
peratures, ie., for T <C T* , on the local- moment side of 
the Kondo phase transition. Experimentally, this would 
require measuring the total susceptibility well below the 
superconducting T c on the undcrdoped side, and has not 
been performed so far. 



2. T matrix 

We now turn to dynamic properties associated with 
the Kondo impurity, in particular we are interested in 
the conduction electron T matrix. The knowledge of the 
T matrix is important for calculating the local density of 
states around the impurity, as observed in STM exper- 
iments. In the Anderson impurity model, the T matrix 
is [Connected to the impurity Green's function according 

tcfea 



T(w) = V 2 G imp (cu) . 



(18) 



We can define the T matrix spectral density as px(w) — 

— \lD.T (bS) I -K . 

The behavior of the impurity spectral function has 
been studied in the particle-hole symmetric pseudogap 
Anderson model in Ref. |4| In this case the strong cou- 
pling behavior is different from the present asymmet- 
ric model, but the asymptotic local-moment behavior is 
identical - this is also seen in Fig. [I], where the particle- 
hole asymmetry flows to zero in the LM phase. For the 



LM regime it was founJi that the impurity spectral 
density vanishes at the Fermi level as |w| r , this could 
be verified in the present calculations. In the asym- 
metric strong-coupling regime we find similar behavior, 
ImGi mp (w) oc pt{uj) oc |cL>| r . Interestingly, the prefactor 
of this power law is equal for both positive and negative 
frequencies, meaning that the spectral asymmetry disap- 
pears in the scaling limit of ImGi mp and px. 

Of particular interest is the behavior of px in the 
quantum-critical region. From general scaling argu- 
ments, we expect a power law px oc lj' ?t ~ 1 , where ??x 
is the anomalous exponent. Our NRG results are con- 
sistent with such a power law, with 77X = 1 — r. So we 
have 



Pt(uj) oc — 
pxM oc 



win (A/w) 



(r < 1), 
(r=l). 



(19) 



Similar to the LM and ASC behavior, any asymmetry 
drops out in the scaling limit of px, and the fixed-point 
spectral density is particle-hole symmetric. 

As done for the susceptibilities, the cross-over between 
the quantum-critical and the LM or ASC behavior in the 
T matrix can be described by a universal scaling function: 



T(uj) = — $ T ( — , — 



(20) 



where A is an amplitude prefactor, and <I>x is a universal 
scaling function (for the particular value of r and for each 
side of the transition). For r = 1 we have to keep in mind 
that logarithmic corrections spoil scaling, and a form like 
( f20| ) is only approximately valid. 

Numerical results for the T = impurity Green's func- 
tion for r = 1 are shown in Fig. ^. They correspond to 
parameter values very close to the transition, and rep- 
resent the scaling functions <&x if we neglect logarithmic 
corrections. For w = u/T* < 1 we have ASC behavior 
with Im<I>x ~ \u>\, for u> ^> 1 the spectral density fol- 
lows the quantum-critical "power law" |a)| _1 ln _2 (A/u)) 
as predicted above, with particle-hole symmetry in both 
limits. However, in the cross-over region the asymme- 
try leads to a large peak for one sign of Q (depending 
on the sign of the overall particle-hole asymmetry of the 
model) - this is true for both sides of the quantum phase 
transition. 

In other words, the Kondo peak known from the metal- 
lic host Kondo model which is located at the Fermi level 
(in the scaling limit) is replaced by a peak at a finite 
energy - this energy actually corresponds to the cross- 
over temperature T* between the quantum-critical and 
the LM or ASC behavior. For the cross-over from critical 
to ASC behavior, one can identify T* with the "Kondo" 
temperature of the problem, since a cross-over from a 
Curie divergence to a constant is seen in Ximp (T) around 
T* . In this paper, we have used the location of the peak 
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FIG. 4. Zero-temperature spectral densities px for the 
conduction electron T-matrix at r = 1, obtained by NRG. Up- 
per panel: J > J c , cross-over from the ASC behavior at low 
energies to quantum criticality at high energies. Lower panel: 
J < J c , i.e., cross-over from the LM to quantum critical be- 
havior. For comparison, the vertical axis has been normalized 
to the maximum value of pt on the "non-peak side". Insets: 
Same on log-log scale, showing both Hi > (solid) and uo < 
(dashed). Note that Q > and Q < are interchanged when 
the overall particle-hole asymmetry changes sign. If we ignore 
logarithmic corrections, the functions shown are equivalent to 
the imaginary part of the scaling function <E>t. 

in pt(^) (Fig. Q) as definition of T* , as it is the numeri- 
cally most convenient and precise criterion. 

We note that very similar results for the T matrix have 
been found in a recent dynamic \axae-N study of a multi- 
channel pseudogap Kondo models. There it was possi- 
ble to analytically determine the low-energy behavior of 
the quantities of interest, both at T = and finite T. 
Related observations have alsOirbeen made within slave- 
boson mean-field calculationsEBEI 

We have also calculated the finite-temperature be- 
havior of the impurity T matrix by NRG, for a fixed 
(temperature-independent) host density of states. It is 
worth emphasizing that the finite-T NRG method is a 
unique tool for this task, as the commonly used slave- 
boson method shows a spurious finite-T transition. A 
series of cross-over functions for the spectral intensity is 
shown in Fig. ^. As expected, the peak is broadened, and 
looses spectral weight for T ^> T* , whereas the peak posi- 
tion is nearly unaffected by changing T. The evolution of 
the weight under the peak (integrated from u> = — 10T* 
to 10 T*) is shown in Fig. similar to the usual Kondo 
effect the weight loss occurs rather slowly, e.g., the weight 
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CO / T* w/T* 
FIG. 5. Same as Fig. ^, but now for finite tempera- 
tures and J > J c . The four panels show the evolution of 
the cross-over peak for temperatures T — T*/2, T* , 2T* , AT* , 
the dashed line is the T = result. The intensity has been 
normalized to the height of the T = peak. The "Kondo" 
peak is broadened and slowly looses weight for temperatures 
above T*. 

is reduced below 50 % of its T = value only for T > 5T* . 

The results in Figs. || and || are possibly of relevance 
for the STM experiments done on Zn-doped BSCCCoeI 
- the large peak in the differential conductance close to 
zero bias correspond to the Kondo peak arising from the 
screening of the Zn- induced moment. We will discuss this 
issue further in Sec. M. 
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FIG. 6. Weight of the cross-over peak as shown in Fig. ^, 
integrated from uj — — 10T* to 10 T* , and normalized to the 
peak weight at T — 0. 
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C. Location of the critical point 

So far we have discussed magnetic screening in the 
pseudogap Kondo model for a fixed host density of states. 
Motivated by NMR experiments on Zn- and Li-doped 
cuprates we now consider the dependence of the pseudo- 
gap Kondo effect on the hole doping level, i.e., with vary- 
ing concentration of mobile carriers in the Cu02 planes. 
Experimentally, a strong dependence of the Kon d o t em- 
perature on the doping level has been observecEaa, in 
particular Tk appears to vanish in strongly underdoped 
samples. 

As explained in the introduction, a (boundary) transi- 
tion between an unscreened moment for small doping and 
a screened (or even absent) moment for large doping can 
be anticipated on theoretical grounds. We propose that 
this transition corresponds to the quantum phase transi- 
tion discussed above in the pseudogap Kondo model (see 
also Ref. ^), and we shall show numerical data indeed 
supporting that the transition in the pseudogap Kondo 
model can be driven by changing the doping level. 

Doping in the cuprates {e.g., adding or removing oxy- 
gen in YCBO) is known to have a number of important 
effects, the main one of course being the change of the 
carrier concentration. For the pseudogap Kondo physics, 
the local density of states is of particular importance. If 
we stick to the simple, doping-independent, tight-binding 
band structure used so far, then doping has two effects: 
(i) the change in the chemical potential, influencing car- 
rier concentration and band asymmetry, (ii) the change 
in the magnitude of the d-wave gap as seen, e.g., in tun- 
neling and photoemission experiments. As we will see, 
point (ii) is most important for the pseudogap Kondo 
physics, since a change in gap size immediately influences 
the low-energy part of the host DOS. .— ■ 

For BSCCO it is established from both ARPESEj and 
tunnelingEj measurements that the magnitude of the su- 
perconducting gap (as well as the normal state pseudo- 
gap) changes considerably with doping: in contrast to T c 
the gap increases in the underdoped regime, there taking 
up to twice of its value at optimal doping. For YBCO 
the situation is less clear since high-quality tunneling or 
ARPES data are rareO. At least, thermodynamic mea- 
surements indicate an increase of the pseudogap tem- 
perature with decreasing doping also in YBCO, there- 
fore a decrease of low-energy spectral weight similar to 
BSCCO is likely. In the following, we will use thii.daping- 
dependent gap values as measured in BSCCCO'Ej, but 
we expect the results to reproduce the correct trend for 
most cuprates. The precise gap values used are A = 57, 
38, 22 meV for doping 5 = 10, 15, 20 %. The host band 
structure entering the Hamiltonian ([l]) is assumed to 
have a doping-independent tight-binding form including 
longer-range hopping terms, we have performed calcula- 
tions with the three parameter sets quoted in Table I of 
Ref. |5l| as well as a simpler dispersion with t = —0.15 eV, 
t' = —t/A, t" — t/12, with qualitatively similar results. 
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FIG. 7. Low-energy part of the input DOS used for the 
NRG calculations, here without potential scattering, (7 = 0. 
Solid lines: 10% underdoped. Dashed lines: 19% overdoped. 
a) DOS for a single-site impurity, being equivalent to the local 
DOS of the d-wave superconductor, b) Effective DOS for a 
four-site impurity as discussed in Sec. M. 

The employed densities of states are shown in Fig. |?]a 
for underdoped and overdoped hosts; Fig. Qb shows the 
effective DOS seen by a four-site impurity, see Sec. 

Of course one has to keep in mind that the evolu- 
tion from the optimally doped ci-wave superconductor 
to an insulator in the underdoped limit involves Mott 
physics which is not included in the simple one-particle 
description employed here. Therefore, the depletion of 
low-energy spectral weight, here modelled by an increas- 
ing ci-wave gap, should be viewed as a phenomelogical 
account for strong correlation effects. 

NRG results for the doping dependence of the criti- 
cal J c are shown in Fig. ^, for both a single-site impurity 
and a spatially distributed four-site impurity (see Sec. ^). 
The data show a 30 to 50% variation of the critical J, 
induced by the size change of the superconducting gap. 
For comparison we have also calculated the critical J 
with doping-independent gap - the effect of the change 
in asymmetry is very small due to the presence of a finite 
asymmetry already at zero doping arising from hopping 
processes beyond nearest neighbors. In addition, we have 
checked that the change in the form of the superconduct- 
ing gap as reportecS for BSCCO, namely the deviation 
from the (cosk x — cosk y ) form for underdoped samples, 
has only a weak effect on the critical J. We note that ad- 
ditional potential scattering terms of course modify the 
J c values, but we have verified that for moderate val- 
ues of the bare scattering potential U the trend that J c 
increases for strong under-doping is preserved. 

The main message of Fig. ||is that J c acquires a signif- 
icant doping dependence, and has a value in the range of 
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FIG. 8. Critical value of the Kondo coupling J c vs. dop- 
ing, extracted from NRG calculations for a cuprate host DOS, 
with a doping-dependent superconducting gap. Solid: sin- 
gle-site impurity. Dashed: four-site impurity (see Sec. [v]). 

magnetic couplings which are expected in the cuprates. 
Of course, the precise form of the microscopic Hamilto- 
nian describing the induced moment is not known, and 
therefore no estimate of the coupling value can be given. 
If we just assume a doping-independent Kondo coupling 
of a certain size, then we can obtain the values of the 
cross-over temperature T* from our NRG calculations, 
and sample data are shown in Fig. |^. The important 
result is that a smooth variation of the input parame- 
ters can lead to a strongly doping-dependent cross-over 
("Kondo") temperature arising from the vicinity to the 
boundary quantum phase transition. Therefore, the ex- 
perimentally observed strong doping dependence of the 
NMR Weiss tcmperatureEj'E3 is not necessarily an indi- 
cation of a drastic change in the host electronic structure. 

We close this section with a remark on the analytical 
understanding of the doping dependence of J c as shown 
in Fig. |[ Certainly, the value of J c will be mainly in- 
fluenced by the low-energy part of the host density of 
states, but one can make this more precise, trying to 
answer the question "which" electrons contribute to the 
Kondo screening the pseudogap Kondo model. A rough 
estimate is provided by the expression for the critical cou- 
pling within the slave-boson mean-field approximation, 

- ex / — po(w) , (21) 

Jc J W 

with Po(lo) cx \uj\ r being the host density of states. This 
integral is logarithmically divergent for r = 0, indicating 
that the moment is always screened as T — > 0, and that 
this screening arises from host electrons close to the Fermi 
level. For r > 0, it is clear that all conduction electrons 
contribute to screening, with a weight proportional to 
1/lu. Therefore, an estimate of J c based on the low- 
energy part of po only may not be appropriate. We note 
that a similar expression forpX, can be obtained within 
a dynamic large- N approach£3 to the pseudogap Kondo 
problem - there, us in the denominator is replaced by 
w 1_Q , where a describes the anomalous exponent of the 
auxiliary fermion propagator at the critical point (with 
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FIG. 9. Doping dependence of T* de- 

scribing the cross-over from quantum-critical to asymmetric 
strong-coupling behavior, calculated with the band structures 
and gap values used for Fig. ^. For the impurity suscep- 
tibility the displayed T* corresponds to the cross-over from 
Curie-like behavior to a constant. Solid: single-site impurity 
with J = 0.14 eV. Dashed: four-site impurity with J = 0.09 
eV. 

a — > for both r — > and r —* 1, but non-zero in 
between) . 

D. Behavior above T c 

The numerical calculations in this paper are restricted 
to the ci-wave superconducting state, but we want to give 
a brief comment on the normal-state impurity behavior , 
The NMR experiments in impurity-doped YBCOEM 
show enhanced Kondo screening above T c - this is rather 
natural as the superconducting gap disappears, the DOS 
near the Fermi level increases, and the metallic Kondo 
effect should be recovered. 

Interestingly, no change in the magnetic properties is 
seen for underdoped samples when the temperature is 
tuned through T c , i.e., the NMR susceptibility follows 
the same Curie law below and above T c . Thus, there is 
no Kondo screening even above T c , which must be in- 
terpreted as a consequence of the celebrated pseudogap. 
The fact that no change at all occurs at T c means that 
the depletion of low-energy spectral weight in both the 
superconducting and the pseudogap phase is similar - 
this appears to support the speculation that pseudogap 
and superconducting gap are of the same origin. 



V. EXTENDED MAGNETIC IMPURITIES WITH 
POTENTIAL SCATTERING 

After having discussed general properties of the pseu- 
dogap Kondo model in the context of cuprate supercon- 
ductors, we turn to a more detailed spatial modelling of 
non-magnetic impurities inducing a magnetic moment. 

One— .j»iportant observation from NMR experi- 
mentsOEJ is that these magnetic moments are not lo- 
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calized at the impurity site, but mainly on the nearest- 
neighbor copper sites. The simplest model for the in- 
duced moment is therefore a Kondo model of the form 
d), 



= ^2 JrS • Sr . 



(22) 



R 



where the moment is coupled to the four Cu sites R ad- 
jacent to the impurity, with Jr = J/4. 

An important point is that such a four-site Kondo 
model (|7|) is not a single-channel model, but instead has 
four screening channels associated with s, p x , p y , and 
rf-wave-like linear combinations of the conduction elec- 
trons on the four sites. These channels are not equivalent, 
and on general grounds one expects that the low-energy 
physics is dominated by the strongest screening channel. 
As our NRG cannot deal with a multi-channel model, we 
have calculated the critical J values for the four channels 
separately, and found that the d-wave channel has the 
lowest critical J or, equivalently, the highest T* for fixed 
J. The physics being dominated by this channel can be 
interpreted as the screening cloud having d-wave symme- 
try Actually, the same conclusion has beetureached in 
a slave-boson treatment of the same modelcj. We note 
that this d-wave symmetry is not directly related to the 
d-wave symmetry of the superconducting order parame- 
ter of the host - it is rather a band structure effect and is 
completely determined by normal host Green's functions. 

We have used such a four-site Kondo model in the d- 
wave channel to calculate the dashed curves shown in 
Figs. U and ^. The effective bath density of states seen 
by the impurity in this case is given by 



P< 



>ff( W ) = -~Im S~] (pRfR'Tl 

7T * 



R.R' 



G(R, R'.w) 



1 



(23) 



with ip R = +[-]l for R-r = (±1, 0)[(0, ±1)], and G is 
the ^ Green's function (||) in Nambu notation. Examples 
are plotted in Fig. - is it seen that the low-energy part 
close to the Fermi level is suppressed (actually, it follows 
an \oj\ 3 power law), whereas the weight at energies around 
the superconducting gap is strongly enhanced. This is 
easily understood if we consider the d-wave-like coupling 
of the impurity in momentum space: the "form factor" 
is (2 + cos 2k x + cos 2k y — 4 cos k x cos k y ), this expression 
vanishes along the diagonals of the Brillouin zone, but 
has maxima at (jr, 0) and (0,7r). 

A second ingredient for an effective model describing 
non-magnetic impurities is a potential scattering term 
(||) . This scattering potential is certainly well localized at 
the Zn or Li site, however, its bare value U is not known. 
Although certain experiments indicate that Zn behaves 
like a strong scatterer, this does not necessarily imply 
a large value of the bare scattering potential U, as the 
discussed Kondo screening can lead to strong scattering 
for low temperatures (note that Tk is much larger in the 



normal state, and can reach 150 K or more). Our results 
so far have been obtained with zero potential scattering 
U; we have checked that the results are practically un- 
changed for U values up to t/4, and remain qualitatively 
valid for U up to 2t. For larger values of the bare po- 
tential scattering, the critical J value for the pseudogap 
Kondo effect is significantly reduced, due to the increas- 
ing density of states near the Fermi level induced by the 
nearby scatterer. 



A. Local density of states and STM 

Motivated by the recent observation of impurity- 
induced |-jf£Sonances in STM experiments on Zn-doped 
BSCCOuu, we can calculate the local conduction elec- 
tron density of states around the impurity. A possi- 
bility first discussed in Ref. [l2] is that the huge peak 
close to zero bias arises from the Kondo screening of 
the Zn-induced moment; related studies can be found in 
Rcf. |3?j. [In the STM context, it is important to distin- 
guish between a screened moment which acts as a spin- 
independent resonant scatterer at low temperatures, and 
an unscreened moment (as in the case of Ni) which leads 
to spin-dependent scattering and causes a splitting of the 
impurity-induced peak in the local DOS.] 

The local (tunneling) DOS in the presence of both a 
potential scatterer and a Kondo impurity is obtained as 



Pstm(R, <*>) = ImTr 

7T 



G(R,R,w) 



1 



(24) 



where G is the full \& Green's function 
G(r,r» = G(r,r» 



+ ]T <fsVs'G(r, s, u;)t z T(cj)t z G(s\ r', u) (25) 



with T(lu) being the Kondo impurity T matrix, and the 
s sites are the neighbors of the impurity at ro. nn 

The STM experiments on Zn-doped BSCCOffil show 
for each impurity a huge peak close to zero bias (at 
roughly —1.5 meV) at a single location which has to be 
identified with the impurity site. This peak signal decays 
rapidly when moving away from the impurity, in addition 
it shows a checkerboard- like spatial modulation, i.e., it 
does not occur on the nearest neighbors of the impurity, 
but on the next-nearest neighbors and so on. This spatial 
shape is hard to explain with the assumption of strong 
potential scattering on the impurity site, because a large 
bare U would expel all electrons there and lead to a very 
small tunneling DOS exactly at the impurity site. 

When comparing the results of microscopic calcula- 
tions with STM data, one has to think about the actual 
tunneling path of electrons between the STM tip and the 
Cu02 plane. BSCCO crystal are cleaved in a way that a 
BiO layer forms the sample surface, and the electrons pre- 
sumably tunnel via this BiO plane into the Cu02 plane of 
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FIG. 10. Calculated tunneling density of states for the 
four-site Kondo impurity model at 15% hole doping with a 
realistic band structure (t = -0.15 eV, t' = -t/4, t" = t/12), 
A = 0.04 eV, and fi = —0.14 eV. The Kondo coupling is 
J — 0.09 eV, the potential scattering U — 0. Top: Local 
DOS vs. energy for the impurity site (red) and the nearest 
(blue) and second (green) neighbor sites. Bottom: Spatial 
dependence of the local DOS at ui = —2 meV. Left: Local 
DOS in the Cu02 plane. Right: Local DOS after applying 
the filter effect proposed in Ref. B5l 



interest. The nature of this tmiiielixig path is not known. 
It has been proposed recentlyoS that tunneling from 
the STM tip into a certain Bi orbital actually probes the 
electrons on the neighboring Cu 3d orbitals, leading to 
a strongly momentum-dependent tunneling matrix ele- 
ment. This so-called filter effect resolves the discrepancy 
between the observed spatial shape on the Zn resonance 
and the result expected from strong potential scattering. 
More precisely, the tunneling matrix elements proposed 
in Ref. |65| have a <i-wave shape, such that the STM signal 
is given by the modified local density of states 



PSTM 



(R,cj) = Im ^r'^r." 

R',R" 

1 +T 2 



Tr 



G(R',R",w)- 



(26) 



replacing (|J), with <^r/ = +[-]l for R' R = 
(±1, 0)[(0, ±1)]. [Note the similarity to the effective DOS 
seen by the four-site impurity, eq. (p3|).] At present, it 
is not clear whether such filtering actually occurs, and 
further experiments are needed to clarify the tunneling 
path into the Cu02 plane. 

The paper Ref. ^ ignored this filter effect (as well as 
the possibility of strong potential scattering at the impu- 
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FIG. 11. Same as Fig. [l^, but with potential scattering 
U = \t\ = 0.15 eV. Here, J = 0.065 eV. The lower panel shows 
the local DOS at u = +2 meV. 



rity site) - we will show here that the results of Ref. [42| 
regarding the spatial shape of the impurity resonance 
are essentially unchanged for weak to moderate potential 
scattering, but change qualitatively for strong potential 
scattering. In the following, we will display the results 
of our calculations bothjuith and without accounting for 
this possible filter effectL^I. 

Fig. 10 shows the calculated tunneling spectra for a 
four-site Kondo impurity with zero (or weak) potential 
scattering. These data are similar to Fig. 1 of Ref. 
but here the impurity properties are calculated using 
NRG removing some artifacts of the slave-boson method, 
namely too large values of the critical coupling J c , too 
large Kondo temperatures for J > J c , and too sharp 
peaks in the T matrix. It is clearly seen that only the 
data in the left panel, i.e., without the filter jeffect, are 
consistent with the experimental observation!], in that 
they show a large peak at the impurity site. Remark- 
ably, the so-called coherence peaks in the tunneling DOS 
at ±Ao are almost completely suppressed near the im- 
purity (cf. Fig. for the bulk DOS) due to the in- 
terference between the bulk Green's functions and the 
impurity T matrix, although the superconducting order 
parameter is not changed (we did not account for gap 
relaxation). Therefore the interpretation of suppressed 
coherence peaks in terms of locally suppressed supercon- 
ductivity has to be used with care. 

Switching on moderate potential scattering on the im- 
purity site does not qualitatively modify the picture, as 
shown in Fig. 11 , Note that for the chosen band structure 
and potential scattering value the global particle-hole 
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FIG. 12. Same as Fig. [j^, but with potential scattering 
U = 4\t\ = 0.6 eV. Here, J = 0.04 eV. The lower panel shows 
the local DOS at u = +3 meV. 



asymmetry has changed its sign, therefore the Kondo 
peak appears at the opposite side of the Fermi level com- 
pared to Fig. [n]. (This should not be seen as a con- 
tradiction to the experiment, as the tight-binding band 
structure is a rough approximation to the real system, 
and the sign of the overall particle-hole asymmetry is 
not known.) 

For huge values of potential scattering, the Kondo ef- 
fect is heavily influenced in the present simple model. 
Due to the large scattering-induced DOS on the neigh- 
boring sites of the impurity the critical Kondo coupling 
is strongly reduced. Correspondingly the Kondo tem- 
perature for reasonable values of J (e.g., 50-100 meV) is 
increased to 100 K or more, which is in disagreement with 
the NMR observation^. Therefore, we consider such a 
parameter combination unlikely to describe the experi- 
mental situation. Nevertheless, we have calculated the 
STM spectrum for such a case, shown in Fig. |l2|. We 
have chosen a smaller value of J, to obtain a Kondo tem- 
perature comparable to the experimental value. It is seen 
that a double peak structure can occur. There are strong 
interence effects between the Kondo peak in the impurity 
T matrix and the potential scattering peak in the host 
DOS, with the result that the local (tunneling) DOS is 
completely dominated by the potential scattering peak. 
The spatial shape of the resulting pattern is more com- 
patible with the experiment after the filter effect is taken 
into account. 

Last not least, in Fig. we show the temperature 
variation of the STM spectra for the situation in Fig. [l(]. 
Besides the simple thermal broadening due to Fermi func- 
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FIG. 13. Temperature evolution of the STM spec- 
trum shown in Fig. |l^, T* w 20 K, calculated 
with finite-temperature NRG under the assumption of a 
T- independent host DOS. Shown is the tunneling DOS at the 
impurity site without applyin g th e filter effect; the dashed line 
is the T = result from Fig. hd. The thermal broadening of 
the signal due to the Fermi functions in both tip and sample 
has been taken into account. 



tions in tip and sample, there is additional broadening 
and weight loss the Kondo peak as shown in Fig. ||. Note 
that we have assumed a temperature-independent host 
DOS, which is approximately justified for temperatures 
up to 0.7 T c . The broadening and the loss of spectral 
weight occur rather gradually, and the peak will be seen 
also for temperature well above Tk- If T* w Tk ~ 20 K 
in optimally doped BSCCO then the broadened Kondo 
peak will survive almost, up to T c . Near T c the gap will 
close, Tk will increaseE3, and the Kondo-induced peak 
might turn into a dip or Faap-like lineshape as seen for 
Kondo impurities in metalsEa. However, this signal will 
be broad and perhaps hard to observe. 

Closing this section, we mention that these-, are sub- 
tleties in interpreting the NMR experimentsU'li-3 in terms 
of a moment located on four sites (which lead to our sim- 
ple effective model). The NMR experiments measure of 
course the fluctuating fields in the full system, whereas 
our model specifies the location of the bare moment. For 
a fair comparison one would have to calculate local NMR 
spectra after having accounted for the interaction of the 
moment with the host fermions (and possibly also with 
the antiferromagnetic spin fluctuations). This is possible 
in principle, but is left for future studies. 
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VI. CONCLUSIONS 

In this paper, we have studied the dynamics of mag- 
netic moments in high-temperature superconductors. We 
have focussed on several aspects of the pseudogap Kondo 
model which describes the interaction of a localized spin 
with fcrmionic quasiparticles which obey a linearly van- 
ishing density of states near the Fermi level. This Kondo 
problem shows a non-trivial boundary quantum phase 
transition as function of the Kondo coupling (or the 
size of the d-wave gap); and we have argued that the 
quantum-critical behavior has consequences for suscepti- 
bility, NMR, and STM measurements in the cuprate su- 
perconductors. In particular, the cross-over ("Kondo'') 
temperature T* , as measured by NMR experiments, can 
acquire a strong doping dependence, due to the deple- 
tion of low-energy spectral weight with decreasing dop- 
ing; and we have proposed that the observed vanishing 
of the Kondo temperature can be explained in terms of 
the quantum phase transition in the pseudogap Kondo 
model. The impurity T matrix which is observable via 
STM shows a large peak at a finite energy which coin- 
cides with the cross-over temperature T* . 

Our calculated STM spectra are consistent with the 
experimental observationta for small to moderate poten- 
tial scattering values, and without taking into account 
the filter effect proposed in Ref. Within our simple 
effective model, larger potential scattering values leadpto 
Kondo temperatures at variance with the NMR resultslij. 
Of course, we cannot exclude that the spatial distribu- 
tion of the Zn-induced moment is different from what 
we have assumed, which would of course strongly affect 
the spatial shape of the STM resonance (but only weakly 
change the magnetic properties). One key experiment to 
decide whether the STM peak is of Kondo origin would 
be to measure its temperature evolution. Here the char- 
acteristic quantity is the spectral weight associated with 
the local Zn resonance. However, the weight loss with in- 
creasing T is rather slow (Figs, g and ^), therefore high- 
quality measurements on impurities with small Tr {i.e., 
small bias of the impurity resonance peak) are required. 

The present NRG calculations for the pseudogap 
Kondo model provide more reliable results than the com- 
mon slave-boson mean-field approximation. This holds 
both for the qualitative behavior of the cross-over scale 
T* near the transition, and the quantitative results for 
the critical Kondo coupling and actual Kondo tempera- 
tures for parameters relevant for cuprates. 

In the present work, we have neglected both the re- 
laxation of the local superconducting order parameter 
around each impurity as well as interactions between the 
impurities. The latter are exponentially small in the im- 
purity density and are therefore unimportant for small 
impurity concentration (in the experimentally accessible 
temperature regime). The gap relaxation will mainly 
lead to a small residual low-energy DOS near the impu- 
rity, which causes Kondo screening at low enough tem- 



peratures even in the underdoped regime. However, the 
induced Kondo temperature is exponentially small in the 
residual DOS, and can be safely neglected. 

Summarizing, our studies support that the NMR ob- 
servations of Curie- Weiss-like behavior associated with Li 
or Zn impurities in high-T c compounds can be described 
by the Kondo screening of the impurity-induced mo- 
ments. The peak in the local DOS in Zn-doped BSCCO 
is likely of the same origin, but a number of experimental 
as well as theoretical issues remain to be resolved. We 
emphasize that the given numbers for Kondo couplings 
and characteristic temperatures can be viewed as rough 
estimates only, since the precise form of a Hamiltonian 
describing the impurity-induced moment is not known. 

Further theoretical studies should address the spatial 
dependence of NMR spectra (e.g., on the Y sites near 
the impurity), the additional coupling to antiferromag- 
netic fluctuations, as well as implications for transport 
measurements. 
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